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We present harmonic generation measurements of the intrinsic nonlinear current density j% of 
YBa2Cu307_i (YBCO) films at temperatures close the T c . Experiments on a range of different 
quality samples allow us to extract the dependence of j2 on the penetration depth of the super- 
conductor. In order to model these results, we calculate the intrinsic nonlinear current response 
of d a .2_ !/ 2-wave superconductors in the Meissner regime in the presence of nonmagnetic impurities 
within the self-consistent T-matrix approximation. 

PACS numbers: 74.25.Nf, 74.20.Rp, 74.72.Bk 



Introduction 

Superconductors with d-wave pairing symmetry should 
exhibit an unusual angle dependent nonlinear Meissner 
effect 0, HQ. This has the effect of significantly alter- 
ing the current density dependence of the superfluid den- 
sity at low temperatures, and dramatically affects the 
nonlinear response of superconductor transmission lines 
or resonators. Recently, intermodulation measurements 
of resonators fabricated from epitaxial YBCO thin films 
revealed a sharp increase in the observed nonhnearity 
at low temperatures 0, H, @i a resu lt that was origi- 
nally predicted by Dahm and Scalapino[3 U. While the 
measured temperature dependence suggests that this ob- 
served nonlinear response has intrinsic origins, the effect 
of impurities on the nonlinear behavior in high T c su- 
perconductors has not been explored systematically by 
experiments. This issue is important since the nonlinear 
response and its unwanted intermodulation products can 
severely limit the applications of superconducting devices 
used for e.g. communication filters. 

In this paper we report on measurements of the non- 
linear current density scale j% determined from third- 
harmonic generation experiments for a range of YBCO 
thin-film samples at a fixed temperature (76 K), but 
which display varying levels of growth-induced defects, 
reflected in the measured penetration depth. The sam- 
ples were fabricated under a variety of growth conditions, 
and showed variations in the penetration depth over the 
range from 0.25 - 0.4 /im at 76 K. Subsequent third- 
harmonic generation measurements on these samples al- 
lowed us to examine the nonlinear current density scale 
j2 as a function of the measured penetration depth, in 
order to explore the effects of impurities on the nonlin- 
ear response. We compare these measured results with 
theoretical predictions for a d x 2 _ y i -wave superconductor 
in the Meissner regime in the presence of nonmagnetic 
impurities. 

For a superconductor, the superfluid density is a func- 
tion of the current density. This manifests itself experi- 
mentally as a dependence on current density of the super- 



conductor magnetic penetration depth. The exact form 
of the current-dependent penetration depth depends on 
the nature of the superconductor, including the details of 
the superconductor energy gap. Since for most cases the 
current-density dependence of the penetration depth A(j) 
is very small, we can use a polynomial approximation for 
this quantity. We expect the change in the penetration 
depth to be an even function of the current density j, so 
the lowest-order term in this low-temperature expansion 
is quadratic: 

(^) 2 = l + (^ 2 ) 2 , (1) 

where Ao is the zero-current penetration depth and ji 
sets the scale of the nonlinear current density (ji is 
sometimes also referred to as the intermodulation crit- 
ical current 9]). In a d-wave superconductor, Eq.Q is 
known to break down at low temperatures due to sin- 
gular contributions to the current from nodal quasipar- 
ticles. The nonlinear current density scale ji character- 
izes the strength of the current-density dependence of the 
penetration depth, and its behavior as a function of e.g. 
temperature and disorder can reveal important details re- 
garding the nature of the superconducting state. It is the 
measured temperature dependence of ]2 extracted from 
nonlinear microwave experiments that reveals the d-wave 
behavior predicted by Dahm and Scalapino0, 0- I n the 
present work we examine the dependence of j2 on the 
penetration depth in order to explore the role of disorder 
in the nonlinear response of high T c superconductors. 

Experiment 

Experimentally, there are a number of different ways to 
measure X(j). However, the large value for j2 that results 
from intrinsic effects 8] means that very large current 
densities are required to produce even very small changes 
in X(j). This means that extrinsic effects due to self- 
heating or vortex motion, for example, may easily mask 
the intrinsic behavior, making the extraction of reliable 
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values for j% extremely challenging. One technique that 
possesses the required sensitivity to small changes in A 
is the mutual inductance technique reported by Claassen 
et q/.|lO|. which used a persistent current trapped in the 
thin-film samples to measure X(j) as a function of dc cur- 
rent. By carefully accounting for heating effects of the dc 
current induced in the measurement coils this technique 
was able to extract an accurate value for j\ . The experi- 
ments demonstrated the expected quadratic dependence 
of penetration depth on dc current density [HJ- 

It has also been realized that a current-density depen- 
dent penetration depth will result in a current-dependent 
inductance per unit length for a superconducting thin- 
film transmission lineQ fill ], or a current-dependent 
film inductance for mutual inductance experiments [f2T| . 
Such a nonlinear inductance will produce mixing or har- 
monic products for single-frequency or two-tone stimuli, 
which can be detected with remarkable sensitivity us- 
ing straightforward spectrum analysis!!?] or lock-in tech- 
niques. Use of such harmonic mixing or harmonic gener- 
ation techniques means that nonlinear effects generated 
by even a very small current-dependent contribution to 
the penetration depth can be detected fairly easily even 
for modest applied current densities. In what follows, j% 
determined from both third harmonic generation exper- 
iments at microwave frequencies, as well as from mutual 
inductance measurements at audio frequencies, is exam- 
ined for YBCO samples that contain differing amounts 
of growth-induced defects, as reflected in the measured 
penetration depth. 

For the microwave results described here, we use pat- 
terned coplanar waveguide transmission lines of different 
lengths and cross-sectional geometries fabricated from 
YBCO films with thicknesses in the range 50-400 nm. 
The YBCO thin films were grown by pulsed laser de- 
position, as well as by reactive co-evaporation. All of 
the films were grown on LaA103 substrates, and all had 
superconducting transition temperatures in the 88-90 K 
range. Some films were intentionally grown under dif- 
ferent growth conditions 01 while others displayed dif- 
ference that were simply the result of random process 
variations. For the microwave transmission line measure- 
ments, we apply a single-tone incident signal (typically 
2-6 GHz), and measure the transmitted signal at the fun- 
damental frequency as well as at the third harmonic fre- 
quency. From the measured third harmonic signal versus 
incident power, along with the transmission line dimen- 
sions and measured penetration depth, we can extract 
i't|15|. For transmission lines of different lengths and 
cross-sectional geometries fabricated from the same thin- 
film sample, we obtained consistent values^] for j 2 . We 
also obtained consistent ji values from results of two-tone 
intermodulation measurements in transmission-line res- 
onators, as well as from analysis of the power-dependent 
impedance derived from resonator measurements per- 
formed as a function of rf power[l5||. The self-consistency 
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FIG. 1: Intrinsic nonlinear current density scale j% versus 
penetration depth A extracted from a range of YBCO samples 
of different quality. Here we show the results obtained from 
the measured third harmonic signals (IP3) at 76 K and the 
mutual inductance (MI) results by Claassen[l2(| at 77 K. 



of the results for j'2 extracted from these many different 
measurements suggests that we are measuring an intrin- 
sic, material-dependent nonlinear effect, as opposed to 
extrinsic nonlinear effects due to e.g. vortex motion or 
heating effects. The results of j'2 measurements at 76-77 
K are shown for a number of different samples in Fig. 

as a function of the measured penetration depth (also 
determined at 76-77 K). 

As mentioned above, mutual inductance measurements 
at audio frequencies can also provide a measure of j'2 |l2| . 
For these measurements, the unpatterned superconduct- 
ing film is sandwiched between two coils driven at a fre- 
quency in the audio range (typically 10 kHz). Such an 
experimental arrangement can be used to derive j% in two 
different ways: (I) by measuring changes in the penetra- 
tion depth as a function of a dc current induced in the su- 
perconducting film, or (2) by measuring the signal gener- 
ated at the third harmonic of the drive signal using a lock- 
in amplifier. The two approaches have been shown [l2j 
to give consistent results for ji at 76 K, and ji mea- 
sured by the former procedure has been shown to agree 
quantitatively with the third-harmonic generation results 
at microwave frequencies I0J. The results of the mutual 
inductance determination of 22 for a variety of samples 
from Rcf. 12| are also plotted in Fig. ^ as a function of 
the measured penetration depth. The agreement demon- 
strated between the microwave-frequency harmonic gen- 
eration measurements and the mutual inductance tech- 
niques for 22 values further supports our assertion that 
the nonlinear effects quantified by j% are intrinsic in na- 
ture. For the films measured by the mutual inductance 
technique, the penetration depth was measured also by 
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mutual induct ance 1 16j. For the films measured by the mi- 
crowave third harmonic technique, the penetration depth 
was measured by mutual inductance and by a microwave 
transmission line technique^]- For both cases, the pen- 
etration depth was measured at sufficiently low signal 
levels to ensure that only linear contributions to the pen- 
etration depth were recorded. In Fig. Q]it is evident that, 
independent of the applied measurement technique, there 
is a clear trend showing that films with smaller penetra- 
tion depths yield greater values for ji- In what follows we 
show that these results can be modelled within conven- 
tional d-wave BCS theory by including the role of defects. 



We have included the suppression of T c and Ao by 
solving Eq. ^1 and JSJ together with the self-consistent 
gap equation jll HilEol Elj 
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d& sin 2 (29) 
2^v/^ + A 2 (6)' 



(7) 



where g denotes the coupling constant and the sum is 
restricted by an upper energy cut-off uj c . 

The superfluid density n s (j s , T) is defined from the 
relationr 
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(8) 



Theory 

The total current density can be written as0,0,0 

duj 
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de 



— jr/2 
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(2) 



where f(u>) is the Fermi function and N±(Q,u>) denote 
the density of states for comoving and countermoving 
quasiparticles: 



N±{e,Lu) = Im 



u> ± Ao(i s /j" c )cos6 



A 2 (e)-p±A (j s /j c )cose] 



(3) 



Here, j s is the superfluid current density and j c the ther- 
modynamic critical current density. In the following we 
will be interested in the current response of dirty super- 
conductors and perform the calculations on the imagi- 
nary axis. Then the associated Matsubara frequencies 
uj n are renormalized and must be determined from the 
following selfconsistency conditions 0: 
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where uj n = (2n + \)ttT denote the bare fermion Matsub- 
ara frequencies. As usual, Y — ni/(TrN ) is the scattering 
rate proportional to the impurity concentration rii, and 
c = cot<$o, where So is the s-wave scattering phase shift 
and Nq is the density of states at the Fermi level. The 
strong scattering limit (unitary limit) is characterized by 
c = 0, whereas for the Born limit c ^> 1. For simplic- 
ity we have assumed a circular Fermi surface and study 
current flow along the nodal direction. The d ;c 2_, / 2-wave 
order parameter relevant for the high-T c cuprate mate- 
rials can be written as 



A(0) = A sin(26). 



(6) 



and in the local approximation, is related to the London 
penetration depth A(j s ,T) by 
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(9) 



where n is the total electron density and Ao denotes the 
zero-current penetration depth. In agreement with Eq. 
(JTJ and © , it is through the expansion of the superfluid 
density n s in terms of j s that ji is defined by the following 
relation 



n s {j s ,T) = n a (T) 



1-1 4- 
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Here, j 2 can be expressed as 



(3(T) 



where 



n,(T) 



2tt 



?j 



de 
,2^ 



A 2 (9) cos 2 6 
(A 2 (9) + £ 2 ) 3/2 



and 



r"K dpi 

p(T) = 2vrA 2 / ^A 2 (9)cos 4 9 

r 27T 



(10) 



(11) 



> (12) 



(13) 



T^Re 



4£ 2 -A 2 (9) 
(A 2 (9) + £ 2 ) 7/2 



In order to model the varying crystal quality of the sam- 
ples used to obtain the results shown in Fig. 2]we vary the 
impurity parameters T, c and use Eq. (|§1)- (|13[) to extract 
the penetration depth A (at zero current) and the nonlin- 
ear current density scale j2- In Fig. [2 we show j^/jc a t 
fixed temperature T versus X(T)/X(T — 0) obtained us- 
ing this procedure. We show the results both for unitary 
scatterers and in the Born limit. We have superimposed 
the experimental data points from Fig. Qby scaling them 
with j c = 150MA/cm 2 and A(T = 0) = 150nm similar to 
the values given in the literature^, 0, 0, H^]. As seen 
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FIG. 3: n s (T) and /3(T) versus temperature T for a range 
of scattering strengths T in the case c = 0. Here, T and 
F are given in units of the critical temperature of the clean 
system T c o and the critical scattering strength P c o at zero 
temperature, respectively. The right panel is similar to Fig. 
6 of Dahm and ScalapinoQ. 



FIG. 2: (Color online) Comparison of the experimental points 
shown in Fig. 0and the theoretical results for the intrinsic 
nonlinear current density scale fa versus the (zero-current) 
penetration depth X(T)/X(T = 0). The four curves corre- 
spond to unitary (c = 0.0) and Born (c = 4.0) scatterers in 
the low (+, x) and high (0,n) temperature limit. 



in Fig. [21 the high temperature results agree reasonably 
well with the experiments. At high temperatures there 
is little difference between unitary and Born scatterers, 
as expected. For the dirtiest samples one can speculate 
whether effects other than the intrinsic nonlinearity be- 
come important. Certainly, the data point at A = 480nm 
(see Fig. ^) seems to fall outside the fit, which may be 
related to the fact that this film was grown on sapphire. 

The nonlinear current density scale 32 has a dramat- 
ically different dependence on A at low temperatures 
T <C T c q, where T c q is the critical temperature of the 
clean system. Provided max(7,T) 3> Aoj s /j c , where 
7 = —lmu>(oj = 0) is the residual quasiparticle width, 
we can discuss this behavior within the quadratic ap- 
proximation of Eq. Q . The functional form of the low T 
results shown in Fig. [21 can be understood in terms of 
Eq. (|llfl from Fig. [31 where we show the temperature de- 
pendence of the superfluid density n s (T) and (3(T) in the 
unitary limit for various T. In the clean limit T = 0.0, 
the /3(T) coefficient diverges at low temperatures due to 
the nodes of the d-wave gap0- With increasing disor- 
der, this divergence is strongly suppressed. At higher 
temperatures, the T dependence of (3(T) becomes much 
weaker. On the other hand, the superfluid density n s 
in Fig. [31 exhibits a nonsingular dependence on T. From 
EcL ljllJI it is clear that the low T divergence of (3(T) even- 
tually will drive ji down. This yields a prediction for low 
T -C T c q measurements valid if the samples with lowest 
penetration depth in Fig. ^are clean enough to exhibit a 
pronounced (3{T) upturn at low temperatures. It is clear 
from Fig. 0that low temperature measurements will pro- 
vide a more stringent test of the intrinsic nature of the 
nonlinear response in these materials, and allow one to 



extract information on the nature of the defects limiting 
the transport. 

In summary, we have presented experimental results 
for the nonlinear current density scale ji extracted from 
a range of different quality YBCO samples. We have 
argued that is consistent with an intrinsic origin. A 
calculation including disorder within the self-consistent 
T-matrix approximation provides reasonable agreement 
with these experiments. 
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